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We propose an analytical model based on diffusion-reaction equation approach for electrochemical 
electron transfer reaction, where the rate is limited by the electron transfer process. The electron 
transfer from an ion in solution to the metal electrode would occur as soon as the energy of the 
orbital on the ion matches the Fermi energy of the metal and a new ion with more positive charge 
is formed. Obviously the ion before elctron transfer and the new ion, which is formed after electron 
has transferred, moves under the influence of different potentials. The coupling between these two 
potentials is assumed to be represented by a Dirac Delta function. The diffusive motion in this 
paper is described by the Smoluchowskii equation. Our solution requires only the knowledge of 
the Laplace transform of the Green’s function for the motion in both the uncoupled potentials. 
Our model is more general than all the earlier models, beacuse we are the first one to consider the 
potential of both the ions explicitely. 


The theory of electron transfer reactions in solution is an important topic [l|. It is conventional to use transition 
state theoretical (TST) approach to the problem. TST provides the rate expression a,s k = where AG is 

the free energy of activation and A is the pre-exponential factor. The TST has the limitation of not taking account 
of the deviation of the reactant population from the equilibrium one For fast reaction (in comparison with 

the rate of attainment of thermal equilibrium ) the deviation from equilibrium population can be very important. 
Homogeneous electron transfer is one of such cases and a diffusion-reaction equation approach is suggested for it [1,0 ■ 
In the electrochemical case, magnitude of reaction rate is very sensitive to overpotential and so it is most likely that 
there are regimes where the deviation from equilibrium population becomes significant. In the following, we investigate 
the problem of electrochemical electron transfer using a diffusion-reaction approach. The case where electron transfer 
is the rate determining step, rather than diffusion from the bulk, is considered here. We use Smoluchowski equation 
to understand the reaction dynamics. We start with an ion, in a polar solvent, near the surface of a metal electrode. 
It has an occupied orbital |a) in it, which may give an electron to the electrode. This orbital interacts with the 
electrons of the metal and as a result of which it is broadened. The energy of the orbital |a) depends on the extent of 
polarization of the surrounding solvent and its value, when the surrounding polarization is in thermal equilibrium as 
assumed to be Ca. Motion of the surrounding solvent will obviously change the energy of the orbital and we denote 
this change by Q. So at any instant of time t, the energy of the orbital may be written as Ca + Q(t). It is fairly easy 
to derive the corresponding Fokker-Planck equation for the probability Pi{Q,t) that the shift has a value Q at the 
time t to be 
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where r^is the time of longitudinal relaxation and A is the solvent re-organization energy. If -|- Q{t) moves above 
the Fermi energy ej’ of the metal, then an electron can be transferred to the electrode. This implies that we are 
considering a metal of a fairly large width in energy and also we assume the rate of electron transfer is independent 
of the energy of the orbital \a), as long as it is above the Fermi level of the metal. Therefore, the probability Pi{Q, t) 
is now given by the following equation (T^ 


dPi{Q,t) 

dt 


2kTX d^PiiQA) 


H-“ K{Q)Pi{Q,t), 

tl oQ 


( 2 ) 


In the above equation has a sink term, K{Q). Qp = ep — ea is that particular value of Q, at which the energy of the 
orbital is equal to the Fermi energy ep. Pi{Q,t) is now the probability that the shift has a value Q at the time t, 
and that the electron has not been transferred to the electrode. So Pi{Q,t) is actually the probability that an ion, 
placed in the reaction zone at the time t = 0, would still survive as the same ion at time t. The Eq.(2) is solved by 
Sebastian et. al. [l^, [I3| only under steady state approximation. In our model we consider the survival probability of 
the new ion (formed after the electron is transferred to the electrode), which is denoted by P 2 {Q,t)- So in our model 
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the relevant equation is 
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In the above equation has a sink term, which we assume to be represented by a Dirac Delta function, K{Q) = 
koS{Q — Qf)- The survival probability Pi{t) is defined using Pi{Q,t) by averaging over the entire range of Q. 


Pi{t) = J dQPiiQ^t). 


Eq.(3) can be re-written as given below 

dPi[Q,t) 


dt 


= CiPi{Q, t) — koS{Q — Qf)S{Q — QF)P2iQ, t) 


(4) 


(5) 


dP2{Q,t) 

dt 


= C2P2{Q, t) - ko6{Q - Qf)Pi (Q, t). 


In the above 
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In the follwoing we provide a general procedure for finding the exact analytical solution of Eq. (5). The Laplace 
transform Vi{Q, s) = Pi{Q, t)e~^*dt obeys 


[s - £i]iPi(Q, s) + ko6{Q - Qf)V 2 {Q, s) = P°(Qo) 
[s - C 2 \V 2 {Q, s) + ko6iQ - QF)PiiQ, s) = 0, 

where = -Pi(QiO) and P 2 {Q, 0 ) = 0 are the initial distributions. 
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koS{Q — Qf) s — C,2 
Using the partition technique [HI, solution of this equation can be expressed as 
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where G((3, s; Qo) is the Green’s function defined by 

G{Q, s- Qo) = {Q\[s-Ci- ko^S[s - £ 2 ]-'^]-'| Qo) (10) 

After simplification 

G{Q, s; Qo) = (Q I [s - £1 - fco'G°(Qf, s; Qf)S]-^ | Qo) , (11) 

where 

G^(Q,s;Qo) = (Q|[s-£2]"'|Qo) (12) 

Now we use the operator identity 

[s - £1 - fco"G°(Qf, s; Qf)S]-^ = [s - £ 1 ]“^ + [s - £i]-^ko^G°(QF, s; Qf)^[s - £2 - ko^G°(QF, s; Qf)S]-^ (13) 

Inserting the resolution of identity I = <iy\y){y\ in the second term of the above equation and integrating, we 

arrive at an equation which is similar to Lippman-Schwinger equation. 


G(Q, s; Qo) = GC(Q, s; Qo) + fco^G5(Q, s; Qf)G 2 {Qf, s; Qf)G{Qf, s; Qo). 


(14) 
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where Gi{Q, s;Qo) = (a; | [s — £i] Qo). We now put Q = Qf in the above equation and solve for G{Qf, s; Qo) to 
get 


G{Q, s; Qo) 
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1 — ko^G^{QF, s; Qf)G^{Qf, s; Qf) 


(15) 


This when substitued back into Eq. (14) gives 


G{Q, s; Qo) 


G?(Q,s;Qo) 


ko^G°{Q, s- Qf)G°{Qf, 5; Qf)GUQf, g; Qo) 
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(16) 


Using this Green’s function in Eq. (9) one can caluclate 7^i(Q,s) explicitely. Here we are interested to know the 
survival probability Pi (t) , which is given by 


Pi{t) 



dQP{Q,t). 


(17) 


It is possible to evaluate Laplace Transform Vi{s) of Pi{t) directly. Vi{s) is defined in terms of P{Q,s) by the 
following equation, 

Vi{s) = (^1 + [1 - fco'G°(Q, s; Qf)G°(Q;^, s; Q^ )] k^G^iQF, s; Qf) ^ dQoG°(Qf , s; Qo)Pi°(Qo)) /{s). (18) 


Erom the above equation we see that 'Pi(s) depends on G 2 (Q_f, s\Qf) which is different from the results of all earlier 
studies The average and long time rate constants can be found from Vi{s). Thus, k^^ = "£ 1 ( 0 ) and kL = —{ 

pole of [1 — fcgG°(Q, s; Q_f)G 2 (Qf, s; Qf)(s)] ^), closest to the origin, on the negative s - axis, and is independent 
of the initial distribution but depends on G 2 (Q_f, s; Qf)- The expression that we have obtained for Pi{s), kj and k^ 
are quite general and are valid for any set of potentials. 
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